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A STUDY OF HYPERSONIC SMALL-DISTURBANCE THEORY 1 

By Milton D. Van Dyke 


SUMMARY 

A systematic study is made of the approximate inviscid 
theory of thin bodies moving at such high supersonic speeds that 
nonlinearity is an essential feature of the equations of flow. 
The first-order small-disturbance equations are derived for three- 
dimensional motions involving shock waves, and estimates are 
obtained for the order of error involved in the approximation. 
The hypersonic similarity rule . of Tsien and Hayes, and 
Hayes' unsteady analogy appear in the course of the develop- 
ment. 

It is shown that the hypersonic theory can be interpreted so 
that it applies also in the range of linearized supersonic flow 
theory. Hence, a single small-disturbance theory, and as- 
sociated similarity rule, apply at all supersonic speeds above 
the transonic zone. 

Several examples are solved according to the small-disturbance 
theory, and compared with the full solutions when available. 
These include flow past a wedge and cone, and determination of 
the initial gradients at the tip of plane and axially symmetric 
ogives. For the axially symmetric ogive it is shown that further 
terms can be found only by using Lighthill’s technique of render- 
ing solutions uniformly valid, and thus the initial curvature of 
the pressure distribution is calculated. It is concluded that on 
a body of revolution described by a power series, the pressure 
distribution and shock wave are also given by power series. 

A brief discussion is given of various additional approxima- 
tions from existing theories. 

INTRODUCTION 

Aerodynamic shapes are ordinarily most efficient when 
they cause the least flow disturbance. For this reason, 
simplified theories based upon the assumption of small dis- 
turbances due to thin bodies have proved to be of practical 
value in analyzing incompressible, subsonic, transonic, and 
supersonic flows. 2 For flows at Mach numbers large com- 
pared with unity, however, the pressure disturbances may 
no longer be small (compared with the static pressure) 
even for thin shapes, so that in tliis sense it has been said 
that no small-disturbance theory exists (ref. 1). However, 
if viscosity can be neglected, the velocity disturbances 
remain small compared with the speed of flight (though not 
compared with the speed of sound), and even the pressure 
changes are small if compared with the dynamic pressure. 
In this sense, therefore, a small-disturbance theory exists, 

1 Supersedes NACA TN 3173, “A Study of Hypersonic Small-Disturbance Theory” by 
Milton D . Van Dyke, 1954. 

2 Throughout, “thin” is used to refer to any body whose streamwise slope is small, and 
so applies to slender fusiform objects as well as flat shapes such as airfoils. 


and the assumption of such small disturbances leads to a 
useful simplification of the equations for compressible flow 
at arbitrarily high Mach numbers. 

Viscosity and heat conduction must be neglected in order 
to have a small-disturbance theory. Otherwise, for example, 
the viscous no-slip condition would introduce velocity dis- 
turbances equal to the speed of flight. In many cases this 
simplification does not destroy the essential features of the 
flow. In other cases, the inviscid theory may serve as a 
basis for including viscosity and heat conduction. Thus, 
recent studies of the hypersonic boundary layer (refs. 2 and 
3), which indicate that viscous effects become essential at 
extreme Mach numbers (say, greater than 15), replace the 
boundary layer by a fictitious solid surface and then utilize 
inviscid theory of the sort considered here. 

At sufficiently high Mach numbers, inviscid flow past any 
given thin object requires nonlinear equations for its de- 
scription. We take this as the definition of hypersonic flow: 
Supersonic flow past a thin body is termed hypersonic if the 
Mach number is so great that nonlinearity becomes an essen- 
tial feature. Thus, the definition of hypersonic flow stands 
on an equal footing with the generally accepted meaning of 
transonic flow at the other extreme of the supersonic range; 
that is, flow at a Mach number so close to unity that non- 
linearity (of a different sort) is an essential feature. These 
two terms — transonic and hypersonic — are most meaningful 
when defined (as here) only for thin shapes. They then 
describe two quite distinct regimes which are, moreover, 
separated by a considerable range of “ordinary supersonic” 
flow in which the transonic and hypersonic nonlinearities are 
unimportant, so that linearized theory can account for all 
significant features of the flow. If one attempts to extend 
the terms to thick bodies, the two separate regimes tend to 
merge, so that one must concede that a flow field can be 
simultaneously transonic and hypersonic. 

It should be noted that the term hypersonic has occasion- 
ally been used in the literature with other meanings than that 
adopted here. Flows are sometimes called hypersonic if the 
free-stream Mach number is simply large compared with 
unity (say, 10 or 5, or even 3) . This is not strictly equivalent 
to the present definition because, in principle, at any given 
Mach number a body can always be chosen so thin that 
nonlinearity is insignificant. However, such extreme thin- 
ness does not arise in practice, so that the two definitions are 
equivalent for practical purposes. Again, Oswatitsch has 
defined hypersonic flow as the limiting condition for a given 
body as the free-stream Mach number tends to infinity 
(ref. 4). This is a limiting case of the present definition and, 
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indeed, Oswatitsch’s similarity rules for thin bodies are 
simply special cases of the more general rules. 

Associated with each of the various small-disturbance 
theories is a similarity rule which connects flows at different 
speeds past affinely related shapes. In the case of linearized 
subsonic and supersonic theory, the similarity rule was fully 
understood only long after the small-disturbance theory was 
in common use. On the other hand, the transonic similarity 
rule was developed concurrently with the small-disturbance 
theory. For hypersonic flows, the simplified theory and 
associated similarity rule were first given by Tsien (ref. 5), 
but were restricted to irrotational flows (and to plane or 
axially symmetric shapes). This is a severe limitation be- 
cause strong curved shock waves and consequent entropy 
gradients and flow rotation are essential features of nearly 
all hypersonic flow problems. This restriction was removed 
by Hayes, who indicated in a brief note (ref. 6) that Tsien ’s 
similarity rule is valid for rotational flows and for general 
three-dimensional shapes. The rule was further extended to 
unsteady motion by Hamaker and Wong (ref. 7). 

As a result of this circuitous development, there exists a 
gap in the hypersonic theory. The similarity rule is known 
for full three-dimensional flows with curved shock waves, 
but the underlying small-disturbance theory has never been 
written down. (To be sure, however, its form is known from 
the analogy with nonsteady flow in one less dimension, which 
was pointed out by Hayes.) For the special case of plane 
flow, and for Mach numbers which are not arbitrarily large, 
this gap has recently been closed by Goldsworthy (ref. 8). 
The published examples of applications of the theory are 
limited to the few special cases which are strictly irrotational 
(e. g., the wedge and cone) or are assumed to be approxi- 
mately so. 

The present paper undertakes a systematic study of the 
small-disturbance theory for hypersonic flow. First, the 
small-disturbance problem is derived by reduction of the 
full equations of motion, boundary conditions, and shock- 
wave relations. The similarity rule and unsteady analogy 
appear in the course of this development. This portion of 
the paper may be regarded as an elucidation of Hayes’ note, 
with estimates obtained for the order of error. Next, it is 
pointed out that to within terms of the order neglected, the 
hypersonic similarity theory can be written in the form of 
the similarity theory for linearized supersonic flow, so that a 
single theory and associated similarity rule cover both 
regimes. Then, a number of special problems are solved 
according to the small-disturbance theory and are compared 
with the full solutions when they exist. 

The symbols used in the text are defined in Appendix A. 

HYPERSONIC SMALL-DISTURBANCE THEORY 

BASIC ASSUMPTIONS 

Consider a three-dimensional body fixed in a steady uni- 
form stream. Viscosity and heat conduction are neglected, 
which implies that shock waves will be approximated by 
abrupt discontinuities. 

The body is assumed to be thin, in the sense that the 
streamwise slope of its surface is everywhere small compared 


with unity. The degree of thinness will be measured by the 
small parameter r which may, for example, be taken to be 
the maximum slope of the body, 3 or its thickness ratio. 
However, for inclined shapes r must be indentified with the 
angle of attack if it is considerably greater than the body 
thickness. 

The free-stream Mach number M is assumed to be so hig h 
that the flow is hypersonic. That is, linearized theory is 
inadequate for predicting the essential features of the flow. 
It is known that linearized theory yields an adequate ap- 
proximation if the maximum body slope is small compared 
with the slop e of th e free-stream Mach cone, that is, if /3r<0 
where /3= y'M 2 — 1. As this ratio approaches unity, linear- 
ized theory grows increasingly inaccurate. Therefore the 
flow is hypersonic if the ratio /Sr is not small compared with 
unity. Since r is small, this means that 13 and, therefore, 
also M will be large in the hypersonic range, so that (3 is 
nearly equal to M. Thus, the criterion for hypersonic small- 
disturbance flow may be written 

T<C1 "I 

, , >- with Afr~l or m 

M»1 ) v ' 

From a mathematical point of view, it is convenient to re- 
gard all the small-disturbance theories as being asymptotic 
forms of the full theory for vanishing by thin bodies. Thus, 
the criterion for hypersonic flow may be expressed more 
formally as : 4 



We introduce a Cartesian coordinate system with the 
positive x axis alined with the free-stream direction (figure 
1). Let the surface of the body be described by B(x, y,z)~ 0, 
and the complete system of shock waves by S(x,y,z)=0, 
where the function S' is not, of course, known at the outset. 



Figure 1. — Notation for hypersonic flow past thin body. 


FULL PROBLEM 

Consider the problem of determining the three velocity 
components u, v, w, pressure p, and density p throughout the 
flow field in the vicinity of the body. The mathematical 
system required is the differential equations of motion 

3 If the slope is high in some small region of the body, as at a slightly blunt leading edge, it 
may be presumed that the small-disturbance theory remains valid except locally. In this 
case t might be taken to be the thickness ratio. 

< Here the order symbols are used in the conventional sense: /(r)=0(l) as t-> 0 means that 
/( r ) remains bounded as r - » 0; and /(r)=0Cg[r]) means /(r)/gf(r) is 0(1); similarly, /( T )=o(l) 
means /(r) vanishes as r— >0. 
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(which govern the flow except at discontinuities) , the Rankine- 
Hugoniot relations across shock discontinuities, and bound- 
ary conditions at the surface of the body and far from the 
body. 

Equations of motion. — The differential equations of mo- 
tion, which express the principles of conservation of mass, 
momentum, and energy, are 


(continuity) (pu) x -\-(pv) v +(pw) x —0 (3a) 

( x momentum) uu x -\-vu v -\- +_p z / p = 0 (3b) 

(y momentum) uv x -\-m v -\-im 1 -\-p v /p= 0 (3c) 

( z momentum) uw x +vw„+ww x +p x /p= 0 (3d) 

(energy) u(plp y ) x +v(p/p y ) v +w(p/p y ) x = 0 (3e) 


(See, e. g., ref. 9, ch. 1) Here subscripts indicate differen- 
tiation, and y is the adiabatic exponent of the gas. The last 
equation actually expresses the fact that the entropy is con- 
stant along streamlines, which for steady flow is equivalent to 
the conservation of energy (ref. 9, pp. 15-16). 

Boundary conditions. — At the body the normal compo- 
nent of velocity must vanish. The unit normal vector at 
the surface is proportional to grad B, so that the condition 
becomes 

q. grad B — 0 

where q is the velocity vector, or 

(tangency) uB x -\-vB„+wB x —0 at B = 0 (4) 

The other boundary condition, which implies that the 
body is flying into still air, may be taken in various equivalent 
forms. For present purposes it is convenient to require that 
all disturbances vanish far ahead of the body: 


(upstream) 


v and w— >0 
P-*Po, 


>- as x —> — <=> 


(5; 


Shock-wave relations.— At a shock wave, conservation 
of tangential momentum leads to the requirement that the 
velocity component tangent to the shock surface be con- 
tinuous. The tangential velocity component (figure 2) is 

given by ~q t = {nX q)Xn, where the unit normal vector n is pro- 
portional to grad S. It is convenient to use brackets to de- 
note the jump in a flow quantity across a shock wave so that, 
for example, [it] is the increase in u through the shock. With 
this notation, the condition of conservation of tangential 
momentum becomes 

(tangential momentum) at S=0 (6a) 


This imposes two independent scalar conditions at the shock 
wave, as physical considerations clearly indicate that it 
should. For plane flow, say, in the x-y plane, the last term 
becomes indeterminate and should, of course, be dropped. 

The remaining shock-wave relations express the con- 
servation of normal momentum, mass, and energy across the 



Figure 2. — Components of velocity at shock wave. 


shock. The magnitude of the velocity component normal 
to the shock is (figure 2) . 


-» _> uS x -{-vS v -\-wSz 

2n=2 ' n — y^+sz+sy 2 


Consequently, the other three shock relations are found to 
be (ref. 9, p. 300), using the jump notation, 


(mass) [p (MiS' I +«/Sy-|-«>&)] = 0 

(normal 

momentum) [p (uS x +vS v +wS x y+ 


(6b) 


(S x 2 + S v 2 +S t 2 ) p]=0 vat S—0 (6c) 
(energy) (uS x +vS u +™S ! y+ 

(sy+sy+sy) £]=oJ (6d) 

So far, these relations are quite symmetrical, remaining 
unchanged if the brackets are taken to denote khe change 
upstream rather than downstream through the shock. A 
definite sense of flow direction is provided only by the 
second law of thermodynamics, which requires that the 
entropy shall not decrease across each shock wave, so that 


(2nd law of thermo.) 



at S=0 


(6e) 


For later use we record Bernoulli’s law 


\ (m 2 +p 2 -Hp 2 )+^ i £=const. (7) 

which holds, with the same constant, throughout the flow 
field. 

FIRST-ORDER PROBLEM 

The full problem is now to be simplified by discarding all 
but leading terms in the body thickness t. This will give a 
first-order hypersonic small-disturbance theory, which can 
be expected to provide a close approximation for thin shapes. 

The reduction is conveniently carried out by introducing 
new independent variables which are of order unity through- 
out the flow field. The form of this transformation is 
suggested by simple examples and limiting cases. For 
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example, the approximate solution for a thin plane wedge at 
very high Mach number has been given several times. 
Pertinent results are that the lateral extent of the flew field 
is some moderate multiple of the wedge thickness, that the 
density never exceeds (y+1)/(t— 1) times its free-stream 
value, and that the relative streamwise velocity disturbance 
and the pressure coefficient are of the order of the square of 
the surface slope. Again, the Newtonian impact theory, 
modified to include effects of centrifugal force, which yields 
the limiting solution for M— °° and 7 = 1 , shows pressure 
coefficients proportional to the square of the thickness. 
The tangency condition suggests that, in 'general, all cross 
wind velocities vary directly with the thickness. Such 
considerations suggest introducing new (barred) independent 
and dependent variables, and redefining the functions 
describing the body and shock-wave surfaces, as follows: 

x=x 

y = \ y ► (8a) 

_ i 

Z=— Z 

T J 


Reduced problem. — This transformation of variables is 
now introduced into the full problem of equations (3) to (6). 
If we discard terms which contain r explicitly, such as r 2 pv; 


in the continuity equation, the differential equations 

become 

(x momentum) 

u x +Tuy+ WU- Z + Pxlp = 0 

(9) 

(continuity) 

Pi+ (po)y+ (p w)l~ 0 

(10a) 

(y momentum) 

Vx+JVy+WVJ+Py/p = 0 

(10b) 

(z momentum) 

W x +VW-y+WWz + Pzl~p = 0 

(10c) 

(energy) (p/p y ) x + v(p/p y )y+w(pfp y )^=0 

(10d) 

I he boundary conditions become 



(tangency) 

B x -\-vBy-\-wB2— 

= 0 at 5 = 0 

(11) 


u, v,w-*0 



(upstream) 

p -> 1/t MV 
P-* 1 

* as x — > — c° 

(12) 


and the jump conditions at the shock waves become 


u=u m [l + T 2 u{x,y,z)] 
v=u„Tv(x,y,z) 


(tangential 

momentum) 


S~ x S 7j S; 


w=u a Tw(x,y,z) 


P=p a ,('YM 2 T 2 )p(x,y,z) 

p = p«,p(x,y,z) 

B=B(jc,y,z) 

S—S(x,y,J) 


(8b) 


(8c) 


The new dependent variables are dimensionless, and the 
new independent variables may also be regarded as dimen- 
sionless if the body is of unit length. 

We tentatively assume that all the new dependent vari- 
ables (u, v, etc.) and the new functions B and S are 0(1) as 
r— *0 for fixed Mr, and that the reciprocals of the new 
independent variables (1 [x, etc.) are likewise 0(1). The 
correctness of this assumption is suggested by examples such 
as those discussed above; its justification will come from the 
consistency of the resulting theory. 

It is important to realize that the notation r=0(l) includes 
the possibility that in the limit v becomes arbitrarily small as 
well as the possibility that it approaches a constant nonzero 
value; only the possibility of its growing arbitrarily large is 
ruled out. For example, the reduced velocity components 
u, v, and w will be identically zero in the region ahead of the 
body. On the other hand, it is definitely implied here that 
in at least some portion of the flow field the reduced quanti- 
ties will not vanish in the limit as t— d). 5 • To be sure, they 
may not all be of order of magnitude unity in the intuitive 
physical sense; for example, for flow past a thin flat wing u 
and v will be moderate multiples of unity, but w will be 
numerically much smaller. 


(mass) [p (S' 5 +aSj-l-'icS'j)] = 0 

(normal _ _ 

momentum) [p(iSj+ r/Sj+w*Sl) 2 + 


(S?+Stfp]=0 


Kat S = 0 


(energy) 


(2nd law 
of thermo.) 


^ {S x -{-vSy+wSi} 2 + 

■ ^ (Sf+S;*) |]=0 



(13) 


The parameters M and r of the full problem enter this 
reduced problem only in the combination Mt, which appears 
only in the upstream condition on p (eq. (12)). In the 
hypersonic range, where 1 /Mr=0(l), the reduced problem 
possesses complete internal consistency. This is most read- 
ily understood by considering first the case where it does not. 

In the range of linearized theory, where 6 Mr< 1, equa- 
tion (12) presents a contradiction. It was assumed at 'the 
outset that p is 0(1) but its value upstream, (1/M 2 r 2 )., is 
then not 0(1). This inconsistency is an automatic warning 
that the reduction breaks down in the range of linearized 
theory. Indeed, it will be seen later that the assumed 
orders of reduced quantities are then actually incorrect, 
and that the reduced equations fail to describe linearized 
supersonic flow. 

In the hypersonic range no such inconsistencies arise. 
Hence, the tentative assumptions regarding orders of re- 
duced quantities are justified a posteriori, as are the sim- 
plifications effected by discarding terms that involve r 2 
explicitly. 


t That is, /(t) =0(t) implies that /(r) is not identically o(t). 


• More precisely, where A/V=o(l) as r-»P. 
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Order of error. — Because terms of order r 2 have been 
omitted in the reduction, first-order quantities will differ 
from their exact values by 0(V) . For the special case of 
plane flow, an analysis similar to the preceding has been 
given by Goldsworthy (ref. 8), who considers only a single 
bow wave: 7 he furthermore confines attention to the range 
Mr~ 1 (which, to be sure, may be the range of most prac- 
tical importance) and finds the error to be 0(1 /M 2 ). This 
is equivalent to the present result in that range. However, 
the present result is more general, holding for arbitrarily 
large values of the similarity parameter (assuming, of course, 
that the assumption of a continuum flow remains valid). 
For example, at infinite free-stream Mach number, the 
error in first-order theory is correctly 0(r 2 ). 

It is interesting to note that the error in the various first- 
order small-disturbance theories decreases progressively from 
0(r 2/3 ) in transonic flow to 0(r) in linearized supersonic flow 
to 0(r 2 ) in hypersonic flow. Therefore, under the plausible 
assumption (confirmed by later examples) that these mathe- 
matical order estimates give a reasonable indication of the 
actual physical magnitude of error, the practical need for 
a second-order solution is seen to be greatest for transonic 
flow and least at hypersonic speeds. 

Unsteady analogy. — A significant feature of the reduced 
problem is that the problem for the streamwise velocity u 
has been uncoupled from that for the other variables. Equa- 
tions (10) to (13) constitute a complete problem for v, w, p, 
and p, which can be solved independent of u. Thereafter, 
u (if required) can be determined from Bernoulli’s law: 

u+l == const. = {y _i )A ,p ? (14) 

Consequently, equation (9) and the first terms of equations 
(12) and (13a) are superfluous and can henceforth be dis- 
regarded. 

As pointed out by Hayes (ref. 6), this first-order problem 
is completely equivalent to a full problem for unsteady flow- 
in one less space dimension. The reduced problem of equa- 
tions (10) to (13) is precisely the full problem of unsteady 
motion in the y-z plane due to a moving solid boundary 
described by H(x,y,5)= 0, where x is interpreted as the time, 
and all other barred variables as the actual physical quan- 
tities. 8 The outline of the moving boundary is given by the 
trace of the original thin shape in a cross-stream plane which 
moves downstream with the free-stream velocity (figure 3). 
For example, the problem of steady hypersonic flow past a 
slender pointed body of revolution is equivalent to the prob- 
lem of unsteady planar motion due to a circular cylinder 
whose radius varies with time, growing from zero at time 
x=0. Hayes has given a physical explanation of this analogy. 

It may be noted that the analogy is similar to that arising 
in the slender-body theory of linearized compressible flow, 
as exemplified by the work of Jones (ref. 10) and Ward 
(ref. 11). There, however, the time-dependent analog is 
incompressible, whereas here it is definitely compressible. 

7 It seems important to consider multiple shocks, because it is not at all obvious that sub- 
sequent shocks have the same status as the bow shock; in particular, at A/=co the upstream 
Mach number is infinite for the bow shock, but has moderate supersonic values for the subse- 
quent shocks'. 

8 The undisturbed fluid has density *><*, =1 and pressurei^ = 1 /yA/V. 



Similarity rule. — The parameters M and r appear in the 
reduced problem only in the combination Mr, which is the 
hypersonic similarity parameter of Tsien. This means that 
for bodies derivable from one another by uniform contraction 
or expansion of all dimensions normal to the stream, the flow 
fields are related if the corresponding Mach numbers are 
such that the similarity parameter Mr is the same in each 
case. The nature of this relationship is simply that the re- 
duced flow- quantities, as functions of the reduced coordinates, 
are the same in each flow. All distinctive flow surfaces, such 
as the body itself and the shock-wave system, have identical 
descriptions in terms of the reduced coordinates. This is 
the hypersonic similarity rule of Tsien, which was extended 
to rotational and three-dimensional flow by Hayes. Its im- 
plications with regard to pressure and force coefficients will 
be summarized later, after the hypersonic rule has been com- 
bined with that for ordinary supersonic flow. 

Extension to unsteady flow. — For simplicity, the preced- 
ing discussion has been restricted to steady flow. It can 
readily be extended to unsteady motions involving small 
time-dependent oscillations of a thin body exposed to a steady 
uniform stream (or, from another point of view, flying through 
still air executing slight time-dependent variations from a 
mean steady rectilinear flight). The full problem is ob- 
tained by replacing a substantial derivative ub/bx-\-vb/by + 
wb/bz wherever it appears in equations (3) to (6) by its un- 
steady counterpart b/bt J t-ub/bx-]-vb/by J rwb/bs. Corre- 
spondingly, if a reduced time t is introduced according to 

1=u^t 

then the unsteady small-disturbance problem is obtained 
from the steady problem of equations (9) to (13) by adding 
b/bi throughout to the operator b/bx+vb/by+wb/bz. For 
example, tlje reduced continuity equation (eq. (10a)) be- 
comes, for unsteady motion 

Pj+Pj-f (p^)j+ (pw);=0 (16) 
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The problem for u remains uncoupled from the problem 
for the other variables. 9 Furthermore, in the remaining 
problem the x and t derivatives appear only in the combina- 
tion d/df+d/c):*;. This means that the number of independ- 
ent variables can be reduced by one by adopting the view- 
point of an observer who is moving with the free stream for, 
to put it another way, who is fixed in still air with the body 
moving past). Thus, introducing “still-air” coordinates 

x'=x — UcJ 

t'=t 

reduces the unsteady small-disturbance problem to exactly 
the form of the steady problem (eqs. (9) to (13)), with t re- 
placed by t' . This means .that Hayes’ analogy remains 
valid for unsteady motion if account is taken of the actual 
variation of the contour with time as well as the apparent 
variation due to relative motion depicted in figure 3. This 
result was first given by Hamaker and Wong (ref. 7). Re- 
cently, Lighthill has analyzed oscillating airfoils at hyper- 
sonic speeds from this point of view (ref. 12). 

If the body oscillates so rapidly or^with such large ampli- 
tude that Bj is much greater than B~, the error remains of 
0(r 2 ) only if r is taken to be the maximum instantaneous 
slope of particle paths at the surface of the body. 

Unified supersonic - hypersonic small - disturbance 
theory. — It would be advantageous if the hypersonic small- 
disturbance theory included linearized supersonic theory as 
a special case. Then the awkward question of what is the 
lower limit of hypersonic flow would not arise. A single 
theory, and corresponding similarity rule, would hold for all 
supersonic speeds above the transonic range. 

In the case of transonic small-disturbance theory, a uni- 
fication of this sort is known to arise quite naturally (ref. 
13, p. 9). In its original form (with /3 2 not replaced by 
2 (M-l)) transonic theory embraces linearized theory as a 
special case, so that it furnishes effectively a unified subsonic- 
transonic-supersonic theory giving a first approximation at 
all speeds below the hypersonic range. 

In the case of hypersonic flow, a connection with the ad- 
joining supersonic range can likewise be effected (ref. 14), 
but the reason therefor is much less straightforward. An 
immediate obstacle is the fact that the approximations lead- 
ing to the hypersonic theory and those leading to linearized 
theory are mutually exclusive. The difficulty arises in the 
continuity equation (eq. (3a)) which in the hypersonic case 
was shown to reduce to 

u„p I J r(pv) v +(pw) z =0 (18a) 

but in linearized theory reduces instead to 

^otP-T - 1 1 pa, (Py 1 W 2 ) 0 (18b) 

The term p m u x must be retained in linearized theory, 
whereas it must be neglected in hypersonic theory in order 
to achieve similitude. It would therefore appear impossible 

9 However, it can no longer be found immediately in terms of the solution of the remaining 
problem, because in unsteady motion there is no useful counterpart of the Bernoulli equation 
(eq. (14)); fortunately, u is seldom actually required. 



to give a small-disturbance theory which is general enough 
to describe both hypersonic and ordinary supersonic flows, 
and yet simplified enough to retain the corresponding simil- 
itudes. It must be regarded as a coincidence that this can, 
nevertheless, be accomplished. 

It is shown in Appendix B that the hypersonic theory 
covers the ordinary supersonic range if it is reinterpreted in 
accordance with the similarity rule for linearized theory 
(ref. 15). It is found that solutions of hypersonic small- 
disturbance theory remain valid at small values of the 
parameter Mr (which is the domain of linearized theory), 
provided that Mr is replaced by j8r and the results are re- 
interpreted in terms of physical variables according to 

u=wjl + r 2 u(x,y,3 ; (3r)] 
v =u^t v(x,y,z; /Sr) 
w—u^t w(x,y,z; /3t) 

V =P=, ^ 1+7MV \j>(x,y,z) = 

Po-{ 7 M 2 r 2 p—J^ 

P=P~^ 1 +^~[V(z>i/-Z; /3r)— lj\ = p ra P— J§) 

rather than according to equation (8). The pressure 
coefficient is given by 

c ^f^r 2T2 (*~ws) (i9b) 

2 Pro^co 

Since the error is 0 (r 2 ) in the hypersonic theory and 0(r//?) 
in the linearized theory, the error in this unified theory is 
0 (t 2 ) or 0 (r//3), whichever is the greater. 10 

The flow quantities of chief aerodynamic interest are the 
pressure coefficient and the various force and moment co- 
efficients derived from it by integration. If the hypersonic 
small-disturbance problem has been solved to find C p /t 2 as a 
function of the reduced coordinates and the parameter Mr, 
the result is rendered valid also in the ordinary supersonic 
range, simply by replacing Mr by /3r. 

Unified similarity rule. — The unified supersonic-hyper- 
sonic similarity rule may be summarized as follows; 

For steady flow past thin bodies derivable from one 
another by uniform contraction or expansion of all dimensions 
normal to the stream, the flow fields are related if the cor- 
responding Mach numbers are such that the similarity 
parameter fir is the same in each case (r being any measure 
of thickness). The relationship is such that the flow fields 
are identical when expressed in terms of the reduced flow 
quantities 

u—u m v w p—p a C p 0 2 ( p 

r 2 u a ’ mj tuJ or t 2 ’ M 2 \ p ro 1 

as functions of the reduced coordinates 


i° For slander shapes such as thin bodies of revolution, the error is only 0(r 2 /£ 2 ) in the lin- 
earized theory and, hence, U(r 2 ) or 0 (t*// 32 ) in the combined theory. 


> (19a) 


AT 



A STUDY OF HYPERSONIC SMALL-DISTURBANCE THEORY 


7 


When the body oscillates slightly, the same is true with the 
addition of u a t to the reduced coordinates, provided that 
the time history of oscillation in terms of reduced coordinates 
is the same for each body. 

From the rule for pressures follows the equality of the 
reduced force and moment coefficients 


Cp C L 

t 2 +*’ T I+*’ 


c , . 

— > etc. 


Here k— 1 if some plan-form area is taken for reference and 
k=0 if some cross-sectional area is used. The connection 
between the similarity rules for forces and moments is con- 
tained in the statement that the center of. pressure is constant 
in terms of the reduced coordinates. 

An unlimited number of equivalent forms of the reduced 
variables can be produced, for example, by multiplying by 
any power of the similarity parameter. However, the forms 
given here are the most useful ones in the hypersonic range, 
because they involve functions of order unity. In the super- 
sonic range no forms have this - advantage, except in the 
special case of plane flow, where it would be convenient to use 

- (— — A and - C R 
r \U„ J t n 

It should be noted that in the hypersonic and combined 
similarity rules the adiabatic exponent 7 must remain 
fixed, whereas its magnitude is arbitrary in the supersonic 
case (since it does not appear in the linearized problem). 


TYPICAL APPLICATIONS OF HYPEBSONIC SMALL- 
DISTURBANCE THEORY 

Several problems will now be solved according to hypersonic 
small-disturbance theory. These examples will illustrate 
possible methods of solution, and demonstrate the degree of 
simplification resulting from the assumption of small dis- 
turbances. Comparison with the corresponding solutions of 
the full equations (when available) will indicate the accuracy 
to be anticipated when the theory is applied to more elabo- 
rate problems. 

PLANE AND AXIALLY SYMMETRIC PLOWS 

The examples to be considered are either plane or axially 
symmetric flows. Accordingly, it is convenient to introduce 
coordinates x,r where in the case of plane flow r is the Car- 
tesian coordinate y. Henceforth, v will denote the velocity 
component in the r direction, which is the radial velocity 
for axially symmetric flows. 

Equations of motion. — In these coordinates, the hyper- 
sonic small-disturbance equations of motion (eqs. ( 10 )) 
become 

Pi+W;+n |-= 0 ( 20 a) 

aj-|-Fy;+-S ^==0 ( 20 b) 

P 



Figure 4. — Notation for plane or axially symmetric body. 


Boundary and shock conditions. — Let the surface of the 
body be described by r=rb( x), as indicated in figure 4. 
Then the tangency condition of equation ( 11 ) becomes 


v=b'(x) at r=b(x ) 


( 21 ) 


The upstream conditions are given by equation ( 12 ). 

The examples to be considered will involve only a single 
bow shock wave, which may be described by r=rs( x). The 
shock-wave conditions are given by equations (13), but in 
this case it is easier to take advantage of the convenient 
relations which are listed in reference 16. Equations (136), 
(128), and (129) of reference 16, when reduced to hypersonic 
small-disturbance form, give 


(Y+l )* 2 
P 2+(T — 1) k 2 


fat r=s(x) 


( 22 ) 


where k=Mts'(x ) is the hypersonic similarity parameter 
based on the local shock-wave slope. 

Stream function.- — The continuity equation (eq. (20a)) 
may be accounted for by introducting a stream function, 
setting 


so that 



.(23) 


- ifc 


•h 

if? 


(24) 


(zVp T );+«( 2 Vp T ) 7=0 ( 20 c) 

The distinction between two and three dimensions arises only 
in the continuity equation, where 0 -= 0 for plane flow and 
v=l for axially symmetric flow. 


Then equation (20c) states that p/p T is a function only of 4*, as 
is clear from the fact that entropy is constant along stream- 
lines between shock waves. Substituting into the momen- 
tum equation ( 20 b) gives 


324350— 65 2 
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'Pr 2 'Pxx — tyzPr'Plf+'Px'pTT = ^ 

where 

“W')=P/p t 

The pressure is given by 

p = cop 7 = W (l/^/r*) "T 
PLANE WEDGE 


^f 7+1 (7coi/'F7+co'i/'? 2 ) for plane flow 

^=r y^+a/Z/J for axially symmetric flow 


(25a) 

(25b) 


(25c) 


(26) 


As a simple introductory example, consider hypersonic 
flow past a thin wedge of semivertex angle 5. Here, and in 
the examples to follow, the solution is most readily carried 
out by assuming a given shock wave and Mach number and 
calculating the corresponding body shape. It is therefore 
convenient to identify the thickness parameter r with the 
shock-wave angle rather than the wedge angle. This is 
quite permissible at hypersonic speeds, where they are of the 
same order; it is also permissible in the ordinary supersonic 
range, where they are not, provided that 8 rather than r is 
used in the error estimates. Let b=8/r be the ratio of wedge 
to shock angles. 

The flow field is conical (velocities constant along rays), so 
that the stream function has the form 


Tp(x,r) =xf(6) 

where 8 is a conical variable, defined by 


® TX x 


(27) 


(28) 


which varies from b at the wedge to unity at the shock wave. 
Thus, equation (25a) becomes 


f"[f-yco 0 f ny+o ]=0 
where, from equations (22) and (25c) 


0)0 = 


2y«o 2 


(7-1) p+(7-l)Ko 2 T 
7(7+1)k 0 2 L (T+1)kq 2 _ 


(29a) 


(29b) 


with k 0 —Mt. The shock-wave conditions of equation (22) 
give (since s' — 1) 


/(1) = 1 
/'(«= 


(Y+1V 


(30) 


2+(7-1)k 0 2 J 

The solution of this problem, which corresponds to/" = 0, is 

_ ( 7 + 1) ko 2 6 — 2 (k 0 2 — 1) 


m= 


2 +(7 — l)«o 2 


(31) 


Requiring / to vanish at the surface gives the ratio of wedge 
to shock angles 

8 2(k 0 2 — 1) 

r (7 + 1) k 0 2 


( 32 ) 


The auxiliary hypersonic parameter k 0 =Mt can now be 
eliminated in favor of MS, with the result that 


1 7+1 
~b~ 4 


-vm 1 


1 


M 2 S 2 


(33) 


Then from equations (26) and (27) the pressure coefficient 
at the wedge (or anywhere between the wedge and shock 
wave) is found to be 


— [pWeF 




(34) 


These results were first given by Linnel (ref. 17). Numerical 
values of pressure coefficient are compared in figure 5 with 
the exact results for wedges of various thicknesses. Here 5 
has been taken to be the tangent of the wedge angle, though 
within the scope of the small-disturbance theory it might 
equally well have been identified with the sine of the angle, 
the angle itself, etc. 



Figure 5. — -Wedge pressure according to hypersonic small-disturbance 

theory. 


Unified supersonic-hypersonic result. — Replacing the 
parameter MS by fio renders the solution valid in the ordinary 
supersonic range as well: 


c v -s 2 3 r+V 0 it) + di 2 ] 


(35) 


This result is again compared with the full solutions in figure 6. 
The advantage of the unified result is obvious. In all sub- 
sequent examples the results will be presented only in this 
form. 
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Figure 6. — Wedge pressure according to unified supersonic-hypersonic 
small-disturbance theory. 

It is interesting to note that this formula has been proposed 
by Ivey and Cline (ref. 18) for predicting the surface pressure 
on any supersonic airfoil. They obtained it by seeking an 
interpolation formula connecting Ackeret’s linearized theory 
with the hypersonic result of Linnell (eq. (34)) for a tangent 
wedge. 

INITIAL GRADIENTS FOR PLANE OGIVE 

Using the full equations, Crocco first determined the initial 
gradients of flow quantities at the tip of a plane ogive by 
perturbing the solution for a wedge (ref. 19). His analysis 
has been repeated and elaborated upon by Schafer (ref. 20) 
and others. This problem provides a good test of the small- 
disturbance theory in a case involving shock-wave curvature. 
Let b be the initial ratio of body slope to shock-wave slope, 
and l the corresponding ratio of radii of curvature. Then in 
physical coordinates the body (figure 7) may be described by 

n— T (bx+lj' cx 2 + . . .^) (36) 

and the shock wave by 


Here c will be negative for the convex shapes usually en- 
countered in practice. 

Shock-wave conditions. — Conditions just behind the 
shock wave are found from equations (22) to be given by 


—4>x=pv 


2W-D f 

2 + (T— 1V\ 


1 


(7-lW+(Y+5W— 2 

V~l) [2 +(y-1K 2 ] 



(38a) 


^ - u+iv r , 4 

W P 2+(7-1)k 0 2 L 1+ 2+(7-1W 


lcx-j- 



(0=p/p y - 


27 k 0 2 -( 7 -i) p+fr-iy r f 
7(7+ IV L (7+l>o 2 J \ 

47(7-1) (no 2 — l ) 2 
[ 27 k 0 2 — ( 7 — 1 )] [ 2 +( 7 - 1 )k 0 2 ] 


1 + 


lcx+ . 



(38c) 


where k q —Mt is the auxiliary hypersonic similarity param- 
eter based upon the initial shock-wave angle. 

Equations of motion. — The flow behind the shock wave 
will consist of a uniform field upon which is superimposed a 
perturbation field duel to body curvature. Hence, along 
each ray from the vertex the flow quantities will have con- 
stant values associated with the initial slope of the body plus 
linear variations proportional to the initial curvature 
(together with higher variations which need not be considered 
in evaluating initial gradients). Thus, the stream function 
of equation (25a) can be written in the form 

t(x,r)=x /(0)— lcx 2 g(d)+ . . . (39) 


Along the shock wave 


i f'=r=x J r . . 

so that the entropy function to can be written as \p(oo) by 
replacing x by i/' in equation (38c). In this form (since u is 
constant along streamlines) the expression applies through- 
out the flow field downstream of the shock wave, so that 


co(t)=co 0 [l+lc<o 1 \p + . . . ] (40a) 


S =T^ x+i lcx 2 + . . 


(37) 



where 


27 ko 2 — ( 7 — 1 ) r 2 +( 7 -lKT 

C °° 7 ( 7 + IV L ( 7 + IV J 

47 ( 7 — 1 ) V— l ) 2 
[ 27 ko 2 — (7-l)][2+(7-lVl 


(40b) 


Substituting these expansions for ^ and to into equation 
(25a) and equating like powers of x yields two ordinary 
differential equations. The first is equation (29a), corre- 
sponding to the basic flow past a wedge. The second, when 
simplified with the aid of equation (31), becomes 

pg" -VJ'g' -uj’ 2 ) (41) 

Solution for g . — By expressing conditions just behind 
the shock wave in terms of conditions at 6=1 through Taylor 
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expansion, the boundary conditions on g are found to be I r T J 


*(i)=- 


K 0 2 — 1 


2+ (7 — 1)k 0 2 


g'(D= 


4(7+1)k 0 2 
[2 + (7 — 1)k 0 2 ] 2 J 


(42) 


The velocities associated with g vary linearly with distance 
from the apex. Hence (as is readily verified from the differ- 
ential equation), g is a quadratic function of 8, and can 
therefore be written as 

5(0)=5(l)-a-%'(l)+^(l-^)V'(l) (43) 


Here only g"(l) is unknown, and it is immediately found 
from the differential equation to be 


£"(!) = 2(7+ l)«o 4 


3(7-1)k 0 4 +(3-7)k 0 2 +2(7+2) 
(ko 2 1) [2 + (7 1) Ko 2 ] 3 


(44) 


Curvature ratio and surface pressure. — The body curva- 
ture and pressure gradient can be expressed in terms of the 
values of g(8) and its derivative at 8=b. From equations 
(42) to (44) it^is found that 


g(by- 


g'(b)- 


2(27 — 1)kq 4 +(7+5)k 0 2 — (7 — 1) 
(7+1)(k 0 2 -1)[2+(7-1)ko 2 ] 

2k 0 2 (3k 0 2 +1) 


(45) 


( Ko 2 -1)[2+(7-1) Ko 2 ; 


The surface of the body is determined by the vanishing of 
the stream function. Thus it is found that the ratio of 
shock-wave curvature to body curvature is given by 


, f'(b) (7+1)V(«o 2 -1) 

l ~2g(b) 2[2(27 — 1)ki) 4 +(7+5)ko 2 (7 1)] 


The initial pressure gradient on the surface of the body is 
given by 


d<7„ 

dx 


—y mf{b) 


. g'(b) 
g(b) 


(47) 


so that in terms of the initial slope R 0 '=br and curvature 
R/ r —C t of the body 


d C P j (7+ 1HW+ 1) [27 *o 2 - (7- 1)] „ , „ „ 

' dx Is (k 0 2 — 1)[2(27 — 1)k„ 4 +(7+5)ko 2 (7 1)] U ° 


(48) 


Although these results have been expressed in terms of the 
auxiliary parameter k 0 =Mt, they can be given explicitly in 
terms of Mach number and apex angle (in the combination 
MRo'=Mbr ) with the aid of equation (33). In this respect 
the small-disturbance solution is superior to the full solu- 
tion, which yields no such explicit results. Replacing 
MR,/ by pR/ renders these results applicable at all super- 
sonic speeds. 

The small-disturbance result for the curvature ratio is 
compared in figure 8 with the full solution for various vertex 
angles, taken from the convenient tabulation of reference 21. 
Figure 9 gives the corresponding comparison for the initial 
pressure gradient. 



P R'o 


Figube 8. — Initial ratio of shook to body curvature for plane ogive. 



Figure 9. — Initial pressure gradient on plane ogive. 

CIRCULAR CONE 

Consider flow past a slender circular cone of semivertex 
angle 8, (figure 10). Again it is convenient to regard the 
Mach number and shock-wave angle as given, and to solve 
for the corresponding cone angle. Thus let the shock-wave 
angle be r and the cone angle o=br, where again & is a 
constant less than unity that is to be determined. 

Equation of motion. — The flow field is conical so that all 
flow quantities are constant along rays. This means that 
the stream function of equation (25b) has the form 

t(x,r) = x 2 f(8), ( 49 ) 
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It follows that the equation of motion (eq. (25b)) becomes the 
nonlinear ordinary differential equation 


The small-disturbance result for surface pressure is com- 
pared with the full results (from ref. 22) in figure 11. 


/' cr-H) / /'\ 

4// , -2// 2 =Tco„ (50) 

Here oj 0 is the constant value of p/p 7 behind the shock 
wave, given by equation (40b), and k 0 —Mt is again the 
auxiliary hypersonic similarity parameter based upon shock- 
wave angle. 

Boundary conditions. — The shock-wave conditions of 
equation (22), together with equation (49), combine to give 

/(D=| 

f (51) 

if / 1 (T+ 1)ktq 

JW 2+(7-1) Ko ^ 


The condition of tangent flow at the cone requires that 
the stream function vanish at the surface: 


/(&) = o 


(52) 


Numerical integration. — The nonlinear equation for /, 
equation (50), can be readily integrated numerically. 
Choosing a value of the auxiliary similarity parameter k 0 , 
we calculate u 0 from equation (40b) and the initial values 
off and f' from equations (51), and then integrate step by 
step inward from the shock wave until / vanishes, which 
determines the cone surface. With the ratio of cone angle 
to shock-wave angle 5/r=6 thus determined, the results 
can be re-expressed in terms of the similarity parameter 
based upon cone angle. Eight or ten intervals between 
shock wave and body yield ample numerical accuracy, 
provided that in each step the predicted values of / and its 
derivatives are corrected by averaging and iterating before 
proceeding to the next step. (For values of k 0 near unity, 
the first few intervals near the shock wave must be taken 
smaller than the others.) 

The pressure coefficient is obtained in terms of the first 
derivative of the stream function according to 



(53) 


where the similarity parameter Mb of the hypersonic prob- 
lem is to be reinterpreted as /35 so that the result is applicable 
throughout the entire supersonic range. 

Computations have been carried out for 7=1.405, in 
order to compare with the full solutions tabulated by 
Kopal (ref. 22). The chosen values of the parameters are 
listed in the following table, together with certain of the 
results: 


Mt 

5/r 

(35 


Cv h /& 

1. 04 

0. 3620 

0. 3765 

1. 217 

3. 183 

1. 19 

. 5545 

. 6599 

1. 522 

2. 646 

1. 58 

. 7281 

1. 150 

2. 207 

2. 333 

2. 87 

. 8604 

2. 469 

3. 925 

2. 154 

4. 47 

. 8921 

3. 988 

4. 973 

2. 116 

OD 

. 9140 

CO 

6. 149 

2. 091 



The differences between the full solution and the small- 
disturbance limit are closely proportional to the square of 
the thickness, in accordance with the estimate of the error as 
0(5 2 ) or 0(<$ 2 //3 2 ) . It is noteworthy that the fractional differ- 
ences are in fact very nearly equal to b 2 . The same is true of 
surface pressures in the previous examples. This suggests 
that the mathematical order estimate may be relied upon to 
give a quantitative prediction of the error in the small- 
disturbance approximation. 

This approximate solution for cones was previously sought 
by Shen (ref. 23), whose result is also shown in figure 11. It 
appears that his solution, which involves more involved 
computations, must contain errors. 

INITIAL GRADIENTS ON OGIVE OF REVOLUTION 

Consider the axially symmetric counterpart of Crocco’s 
problem: the determination of the initial gradients of flow 
quantities at the tip of an ogival body of revolution. This 
problem has recently been considered for the full equations 
by Cabannes (ref. 24) and by Shen and Lin (ref. 25) . It will 
be seen that the small-disturbance solution serves as a useful 
guide to the full solution. In particular, it clarifies the 
behavior of the solution near the surface, which was glossed 
over by Cabannes and in Shen and Lin’s work was mis- 
construed to indicate a singularity which implied wrongly 
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that the initial pressure gradient at the tip of an ogive is 
infinite. 

As in the plane problem; let the body be described in 
physical coordinates by 

r b =T (^> x +\ ca; 2 + • • • ^ (54) 

and the shock wave by 

r,= r (x+^lcx 2 . . (55) 

(figure 7). Conditions just behind the shock wave are given 
by equation (38), with — 1/7 and /y replaced by — 1/7 /r and 

'hi*. 

Equations of motion.— Again along each ray the flow 
quantities have constant values corresponding to the initial 
slope plus linear variations proportional to the initial curva- 
ture (together with higher variations which are considered 
later). Hence, the stream function may be written in the 
form 

t(x,r) = x 2 f(6)—cx 3 g(6)+ . . . (56) 

Along the shock wave 

• • • 

so that an expression for the entropy function w throughout 
the flow field downstream of the shock wave is obtained by 
replacing x with in equation (38c) . Hence, 

«(^0 = . -] (57) 

where w 0 and a>i are given by equation (40b). 

Substituting these expansions into equation (25b) and 
equating like powers of x yields for / the nonlinear ordinary 
differential equation already treated in the problem of the 
cone (eq. (50)), and for g the linear equation 


Behavior of solution near body surface. — Just as condi- 
tions at the curved shock wave have been related to those at 
0 =1, so with the present coordinate system it is necessary 
to relate conditions at the surface of the body to those on the 
initially tangent cone 0=6. However, the solution is non- 
analytic near the surface, which means that Taylor seiies 
expansions do not exist. It is therefore necessary to examine 
the nature of the solution in the vicinity of the body. 

The function f associated with the basic conical flow is 
analytic near the surface (and vanishes at 0=6), so that it 
has a Taylor series expansion: 

/(0) = (0 — 6) /(6)+O[(0— 6) 2 ] (60) 

It follows from equations (58b) and (60) that near the surface 
the coefficients of the differential equation for g behave like 

A~(0-6)- 1/2 ' 

B ~(0 — 6) 

(61) 

C~ 1 

Z)~l „ 

Therefore the point 0=6 is a singular point of equation (58a), 
but is an ordinary (nonsingular) point of the homogeneous 
equation obtained by deleting A (see ref. 26, p. 73). There- 
fore the general solution of the homogeneous equation is 
analytic and can be taken to have the form 

gn=1han(e-by (62) 

0 

where all higher coefficients a r _ can be expressed in terms of 
the two arbitrary constants a 0 and ai by means of the differ- 
ential equation. Then the procedure for calculating a 
particular integral (ref. 26, p. 122) shows that the nonhomo- 
geneous equation (58a) has a particular integral of the form 


Dg"=A+Bg+Cg' 


whose coefficients depend upon / according to 


(58a) 


^=(*-*) 3/2 ZX(0-6)" (63) 

0 


/'( 7+1) 

D=y Uo J 1 ? 7=w -4f 


(58b) 


As in the plane problem, initial conditions on g are found 
by expressing conditions just behind the shock wave in 
terms of conditions at 8— 1 through Taylor series expansions: 


0 ( 1 ) 


Kg— 1 

2+(y — 1 )k 0 2 


► 


0'(1)=-KO 2 


(7-1)ko 2 +2(27+3) 

[ 2 +( Y - 1 ) k „ 2 ] 2 


(59) 


where the coefficients c n can all be determined. Here the 
3/2-power branch point arises from the fact that the pencil 
of fluid striking the tip of the ogive is spread thin over the 
entire surface, and the linear entropy gradient at the tip 
due to a curved shock wave is thus intensified to a square-root 
gradient normal to t ho surface elsewhere. The complete 
solution of equation (58) is the sum of g h and g p . 

In treating the full problem in reference 25, Slien and Lin 
claimed to have found a logarithmic singularity at 0=6, 
which considerably complicated their analysis. Because of 
this singularity, their solution was restricted to concave 
bodies (although they conjectured that it might be extended 
to convex bodies). The singulaj-ity also led to the conclusion 
that for an analytic body shape the shock wave is nonana- 
lytic, and vice versa. Furthermore, the singularity would 
imply that the initial pressure gradient on an analytic body 
of revolution is infinite, although numerical solutions by the 
method of characteristics give no indication of this. 
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The present solution shows no such singularity. It seems 
unlikely that a singularity could have disappeared as a result 
of making the small-disturbance approximation, since this 
would imply that the approximate model does not retain the 
essential features of the full problem. The alternative con- 
clusion is that the singularity does not actually exist in 
the full problem. This has subsequently been confirmed by 
the authors of reference 25" who find that the singularity is, 
in fact, only apparent in the sense of reference 26, page 406. 

It may be noted that Cabannes, in treating the full 
problem, has completely ignored the nonanalytic nature of 
the solution, and simply extrapolated his numerical solution 
to the surface of the body (ref. 24). His results may not be 
seriously in error, because the effect, of the nonanalyticity 
is small. 

Numerical integration. — The differential equation for g 
has been integrated numerically for the six values of the 
similarity parameter chosen previously for the cone. The 
integration was carried out step by step starting from the 
known values at the shock wave and using the same intervals 
as for the cone. This step-by-step solution was joined at the 
two points nearest the surface with the series expansion 
about 8=b given by equations (62) and (63): 


£(0)=~ /W' 2 (0-&) 3 ' 2 [l+| 6-(f? X-±) 

g(b) [l-(-2Xe 3 — X(l + 2X)e 4 + . . . ] + 


* 2 + • • • + 


whc 




( 1+ 2 e 3 X£+ 6 Xe +-") 

1 (64a) 


(64b) 


(Values from the. step-by-step integration and the series 
expansion were also compared at the third point from the 
surface as a check.) Because they are based upon the 
previous solution for a cone, the computations were carried 
out with 7= 1.405. 

Curvature ratio and pressure gradient.— Because the 
nonanalyticity appears only in higher terms of the series, 
surface values of g and its first derivative (but no higher 
derivatives) can be expressed in terms of values at 0=6. 
The surface of the body is determined by the vanishing of 
the stream function. Thus the ratio of shock-wave curva- 
ture to body curvature is found to be given by 


2? (6) 


. (65) 


Proceeding as in the plane case, it is found that the initial 
pressure gradient is expressed in terms of the initial slope 
R 0 ' = br and curvature R 0 "=cr of the body by 


f (by 

dz ~M 2 R 0 ' 2 b«-» g(b) 


(66) 


Numerical values of g(b) and g'(b) are listed in the follow- 
ing table, together with the resulting values for the curvature 


ratio l and surface pressure gradient: 


gRg 

gib) 

g’W 

l 

1 £>CJ 
Ro R o' ' |(, 

0. 3765 

5. 170 

— 8. 502 

0. 0426 

5. 532 

. 6599 

1. 800 

-3. 640 

. 2346 

4. 662 

1. 150 

1. 573 

-4. 073 

. 5106 

4. 514 

2. 469 

2. 101 

-6. 820 

. 8039 

4. 684 

3. 988 

2. 487 

— 8. 638 

. 8921 

4. 802 

oo 

2. 931 

-10. 76 

. 9586 

4. 929 


The curvature ratio is plotted in figure 12, and the initial 
pressure gradient in figure 13. The curvature ratios calcu- 
lated from the full equations by Shen and Lin are also shown 
in figure 12 for comparison, because the error introduced by 
incorrect treatment of the solution near the body is probably 
small. 



Figure 12. — Initial ratio of shook to body curvature for ogive of 

revolution. 

HIGHER TERMS IN SERIES FOR OGIVE OF REVOLUTION 

It might be supposed that for an ogive of revolution de- 
scribed by a power series, the perturbation scheme could be 
continued indefinitely to find successively higher terms in 
a power series expansion for surface pressure. However, 
because of the nonanalytic nature of the stream function 
near the body, complications arise if one proceeds simply 
by adding further terms to equation (56). 

It has been seen that g(8 ) involves an authentic 3/2-power 
branch point near the surface. However, it can be shown 
that the next term will involve a spurious 1 /2-power branch 
point there, the next an inverse 1/2-power singularity, and 
so on. As a consequence of this spurious reinforcement of 
the actual nonanalyticity, it is impossible to evaluate surface 
pressures. Hence, straightforward continuation of the per- 
turbation procedure breaks down. 

The difficulty arises from the fact that in the first perturba- 
tion, the 3/2-power branch point in g arises at the basic cone 
(0=6) rather than at its actual location on the ogive surface. 


ulna private communication; see also Addendum No. 1 to ref. 25. 
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flR'o 


“Wd — “o{ 14'c<>ifc-v^\^+[coi(7nrf+7ic 2 )+co2i 2 c 2 ]^'+ . . .} (70a) 

where u 0 and <oi are given by equation (40b), and 

“* = fo^-l) [2 -Kt-IW] Wl (70b) 

The stream function has the form 

i—xj(d)—cx 3 g(d)—x i [c 2 j(d)+dh(e)]+ . . . (71) 


where 0 is the nearly conical coordinate ?/ x . The differential 
equations which result from substituting this series into the 
equation of motion are simplified by setting 

gCe)=ig(e)~.f' (5) 1 




r ( 72 ) 


jG) =l 2 j(e) +nh(e) +1 k'G) -\f'G) J 


Figure 13. — Initial pressure gradient on ogive of revolution. 

Although this discrepancy has no effect upon the first per- 
turbation, it is compounded in subsequent terms so as to 
be catastrophic. The remedy is to choose a slightly strained 
coordinate system such that for each term the 3/2-power 
branch point appears precisely at the body surface. Both 
the difficulty and the remedy are just those considered by 
Lighthill in his discussion of a technique for rendering ap- 
proximate solutions uniformly valid (refs. 27 and 28). As 
in the previous examples, the solution proceeds, in effect, by 
assuming a given shock wave and determining the corre- 
sponding body shape. Therefore, the required straining of 
coordinates is not known at the outset, but must be deter- 
mined to successively higher accuracy as the solution 
progresses; this is characteristic of Lighthill ’s technique. 

With this modification, the perturbation procedure can, 
in principle, be continued indefinitely. It can therefore be 
concluded that an analytic body of revolution is accompanied 
by an analytic attached bow shock wave at supersonic speeds, 
and an analytic pressure distribution. 

Uniformly valid equations of motion. — Let the body be 
given in physical coordinates by 

r„=T ca;2 +^ d3 ?+ • • (67) 

and the corresponding shock wave by 

r s =r lex 2 - fi (md+nc 2 ) a^-f . . .J (68) 

Now introduce a slightly strained radial coordinate r such 
that the body surface is given by r=bx. The simplest 
choice is 

i— r — | cx 2 — i dx 3j r . . . (69) 


(The functions g, h, and j thus introduced are those which 
would appear in place of g, h, and j if no straining of the 
coordinate system had been undertaken.) Then the differ- 
ential equa tions (and boundary conditions) for j and g, as 
functions of the strained variable 6, are (as implied by the 
common notation) found to be just those solved in the 
previous sections, where / and g were functions of the 
unstrained variable 8. The differential equations for the 
new functions h and j are 


Dh"=E+Fh+Gh' 

Dj"=I+Fj+Gj' 


(73a) 


where D is given by equation (58b), and 


E= 


yUT+n 




r[v(/"-9+y"" 


0 ( 7 - 

F=4(f n +4tf”) 

<?=tw 

[«/ + ^(^')](f-r)- 




(£-/')+/y [^4(^]}]- 3 

2gg'f-2(i' 2 -2gg")f] 


V (73b) 


The procedure which led to equation (57) gives for the 
entropy function behind the shock wave 


Boundary conditions. — Again, expressing conditions just 
behind the shock wave in terms of those at 5=1 by Taylor 
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series expansion and simplifying with the aid of the differ- 
ential equations gives 

H D-gsd) 1 


id)— 


>■ (74a) 




2 (7— l)/tp 2 

[2+(7-1)ko 2 ] 3 J 


where /'(l) is given by equation (51b), gr(l) and g'( 1) by 
equation (59), and 


/"(D 


(Y-1)k 0 2 

2+(7-1)k 0 2 




2k 0 2 [27k 0 2 — ( y ~ 1)1 [2+ (37+ 1)k 0 2 ] 
(7+1)(k 0 2 — 1) [2+(7-1)k 0 2 ] 2 


y (74b) 


[4(7-1) (27+3) + 4(27 3 + 37 2 +5) k 0 2 ~ 

3(7-l)7 2 +(67+l) ko 4 +2(7 — 1) (57+3) k 0 9 ]X 
{(k„ 2 -1) [2+(7-1)k 0 2 ] 3 }- 1 


Behavior near body surface.— The fact that /(0) is 
analytic at 6 = b implies that the full solution for h is also 
analytic, and that the solution of the homogeneous equation 
for j is analytic. The coefficient I is proportional to 
(8 — b)~ 3>2 , so that a particular integral for j is proportional 
to (0— £>) 1/2 . For purposes of computation it is important to 
separate the regular part of I from the singular part, because 
either may predominate for the closest practical approach 
to the singularity, depending upon the value of the similarity 
parameter. Furthermore, the accuracy of joining the step- 
by-step solution for j to the series is increased by treating 
not j itself but the combination j+g'/2l. Hence, the series 
employed is 


outlined previously, with 7=1.405. Because h represents 
actually (like g ) only a first perturbation of the basic conical 
flow 12 and is furthermore regular near the body, it is readily 
determined with ample accuracy. On the other hand, in 
the equation for j, the coefficient / is so strongly singular 
that it was found necessary to replace simple step-by-step 
integration by the more laborious five-term procedure of 
Milne (ref. 29, p. 142). „ The coefficients and boundary 
values are also considerably more difficult to calculate, so 
the integration of j has been limited to three values of the 
similarity parameter, whereas h has been found for four 
values. 

The accuracy of the solution for j suffers from the facts 
that it depends upon the accuracy of the preceding solution 
for g, that one coefficient in the differential equation is 
strongly singular, and that the results of physical interest 
are found as differences of nearly equal quantities. Con- 
sequently, although results derived from the functions /, 
g, and h are probably reliable to three or four significant 
figures, those derived from j arc perhaps not reliable to more 
than two. 

Body shape and pressure distribution. — The parameters 
m and n, which relate the shape of the shock wave to that of 
the body, are found, by requiring the stream function to 
vanish at the surface, to be given by 



The surface pressure coefficient is given by 


^ = _2_ 
7 2 7k 0 2 




(77) 


Thus it is found that on the surface of a body described in 
actual coordinates by 


r=E 0 , x+|f?o /, x 2 +|Bo ,,, x 3 +. ■ • (78) 


(5-j)"’[i+<4( <s+31, V a ) , '‘''’-- • ■}- 

j(b) |^l + 2Xe 2 +-g- Ae 3 +. . .J+ 

j'(b)Q)-b) [l+| X< 2 +^ X‘ 3 +- • •] (75a) 


where /&=/(&), etc., X and e are defined by equation (64b), 
and 


*=b 


(-ftO 


gb 

_2V2* 3 (6+r 2 


37 


gb 


(75b) 


Numerical integration. — The differential equations for 
h and j have been integrated numerically in the manner 


the pressure coefficient is given by 


where 


(79a) 


c, 

R 


7, o 1 1 

o ,2 ~^L& 2 w t(M£„') 2 J 




Ya>o /ft 
b(R 0 'R 0 "x) b \b 


b 2 C P 




2 7C0Q //AY-iA 
d(R 0 'x) b(R 0 "'x) 3 b \b/\ j h ) 


d 2 C D 


d(R 0 "x) 2 


- 27 coq 


mx 


7 ^ o72^> 


4 V 2 6 2 2n t 212 fb 


DJ 


(79b) 


13 It corresponds to a pointed body with zero initial curvature; compare reference 24, 
part III. 
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Numerical values of h(b), h'(b),j(b), and j'(jb) are listed in 
the following table, together with the resulting values of the 
parameters m and n which relate the shape of the shock wave 
to that of the body, and the two functions which give the 
initial curvature of the surface pressure distribution. The 
two pressure-curvature functions are plotted in figures 14 
and 15. (Curves have been faired through the three calcu- 
lated points by analogy with the results of the cone-expan- 
sion approximation discussed in the following section.) 13 


pRv' 

m 

h\b) 

m 

m 

m 


d 2 C p 

&C P 

i 

d(Ro"x)* 

SWx) d(Ro"’x) 

0. 3765 
.6599 

14.34 

1.752 

-36. 40 
-l. 401 



0. 0051 



5. 694 

-7. 34 

43.7 

.0803 

0. 583 

4.04 

4. 736 

1. 150 

. 9341 

-3. 665 

I. 94 

5.30 

.2867 

.868 

5. 20 

4. 561 

CO 

1. 166 

-6. 948 

7.63 

-27 .4 

.8031 

-.901 

13.7 

5.332 
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Figure 14. — First term in initial pressure curvature on ogive of 

revolution. 


FURTHER APPROXIMATIONS 

The theory discussed heretofore is the simplest which 
retains all the essential features of hypersonic flow, so that 
its solutions approach exactness as the thickness tends 
toward zero. Further approximations, although desirable 
for facilitating solution, will introduce errors whose nature 
may be more obscure. In the case of plane flow, however, 
further approximations exist which are so simple and accurate 
that the problem may be considered solved for practical pur- 
poses (cases (5) and (6) below). These and other approxima- 
tions will now be considered for three-dimensional shapes, in 
comparison with the solutions already given. 

Most of these approximations are useful outside the limits 
of hypersonic small-disturbance theory. However, we shall 

13 This three-term series approximation was used to calculate the pressures over circular 
arc ogives for a similarity parameter of 1. The results were compared (ref. 30) with the 
method of characteristics solutions (ref. 31) and shown to be in good agreement over most of 
the ogive. (Note that in Figure 7 of reference 30 the ordinate should be labelled C p /R o' 2 
rather than C P .) 



Figure 15.— -Second terra in initial pressure curvature on ogive of 

revolution. 


consider them here only as they are reduced to small-disturb- 
ance form, so that they actually represent approximations 
beyond those already made. For example, the well-known 
shock-expansion method will be considered only in its hyper- 
sonic small-disturbance form (ref. 1 7) . 

The following additional approximations will be considered: 

(1) Linearized theory, second-order theory, etc. 

(2) Newtonian impact theory 

(3) Newtonian theory plus centrifugal forces 

(4) 7=1 

(5) Cone-expansion approximation 

(6) Tangent-cone approximation 

(7) Compression-layer approximation 

Linearized theory, etc. — The breakdown of linearized 
theory serves almost as a definition of hypersonic flow. 
Hence, the most that can be expected of linearized theory, 
second-order theory, etc., is that they penetrate somewhat 
into the lower end of the hypersonic range. 

For plane flow Donov (ref. 32, pp. 90-91) has determined 
the fourth-order solution. Reduced to hypersonic small- 
disturbance form, his result for surface pressure coefficient 
on a single airfoil may be written as 


M 2 C p =2K+ y -^ K 2 + 7 -- 6 - g 3 + (T+1) - J T 5) K 0 3 + 


3 + 27— 7 2 Ki 3 — 7— 37 2 +7 3 
48 A + 32 




3-77— 77 2 +37 3 
96 


K 0 *K- 


(Y+l) 2 (37-5) 
128 


K 0 3 K 0 'x+ 


(80) 


where K is the local similarity parameter (M times local sur- 
face slope), K n its value at the leading edge, and K 0 ' its 
initial rate of change. Even in this reduced form nothing is 
known of the range of convergence of the series or, indeed, 
whether it converges at all. However, for a single wedge the 
solution is known in closed form from equation (34) . Hence, 


J 
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it is seen that in this special case the series is convergent for 
2f=M5<4/(y+l), which is 1.67 for air. 

For cones the linearized and second-order solutions (ref. 
33 , p. 11 ), as reduced to hypersonic small-disturbance form, 
are shown in figure 16. 



Figure 16.— Further approximations to hypersonic small-disturbance 
theory; pressure on cone, 7 = 1.405. 

Newtonian impact theory. — Assuming that fluid particles 
lose their normal momentum on impact with the surface 
leads to a prediction of pressures proportional to the square 
of the sine of the angle of inclination or, in the small-disturb- 
ance approximation 

§5=2 (81) 

wherever the slope is positive, and zero elsewhere. Accord- 
ing to equation (34) the actual value for a wedge falls only to 
2.4 at infinite Mach number (with 7 = 7 / 5 ), so that the ap- 
proximation is poor for plane flow. It is more satisfactory 
for fusiform shapes such as a cone (figure 16), for which the 
actual value at M= °° (with 7 = 7 / 5 ) is 2.09. 

Newtonian plus centrifugal forces.— Newtonian impact 
theory has been improved by including the centrifugal pres- 
sure gradient through the layer of fluid streaming over the 
body (refs. 34 and 35). The result is precisely the limit of 
the full theory as M—> <» and 7 — >1. In the small-disturbance 
approximation (ref. 33), it gives for plane flow 

C P =2(R' 2 +RR"} (82) 

and for axially symmetric flow 

C P =2R' 2 +RR" 


In both cases it is to be understood that negative values are 
to be replaced by zero. Figure 17 shows that the improve- 
ment due to including centrifugal effects is appreciable for 
the initial pressure gradient on an ogive of revolution. 



Figure 17. — Further approximations to hypersonic small-disturbance 

theory; initial pressure gradient on ogive of revolution, 7 = 1.405. 

7 = 1 . — It has just been seen that on fusiform shapes near 
M — co , the surface pressure is insensitive to the value of 7 . 
At the other end of the hypersonic range, linearized theory 
is independent of 7 . These two extremes suggest that a close 
approximation throughout the hypersonic range may be 
found by setting 7=1 (and this is particularly true since in 
a real gas 7 approaches 1 at high temperatures). This choice 
simplifies the theory by rendering it effectively isentropic; 
that is, ' although shock waves produce entropy jumps, 
entropy does not appear in the pressure-density relation and 
is therefore absent from the problem. 

This approximation has been tested by computing tKe 
hypersonic small-disturbance solution for a cone with 7 = 1 . 
The results corresponding to those tabulated on page 1 1 are 
shown in the following table (including the known value 
at M= 00 ). 


Mr 

5/r 

_ 

PS 

f'(b)/b 

CvJ? 

1. 04 

0. 3912 

0. 4069 

1. 255 

3. 076 

1. 15 

. 5609 

. 6450 

1. 546 

2. 627 

1. 50 

. 7647 

1. 147 

2. 498 

2. 277 

CO 

1 

CO 

03 

2 


(83) 
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Figure 16 shows close agreement with the results for y= 1.405, 
the discrepancy being, indeed, less than that due to the 
thickness of a 10° semivertex angle (cf. figure 11). 

“Cone-expansion” approximation. — The shock-expan- 
sion method for plane flow, which neglects disturbances 
reflected from the bow wave, has recently been shown to 
yield good accuracy at all supersonic speeds away from the 
transonic zone (ref. 36) . 

A more surprising discovery is that an analogous procedure 
yields a reasonable approximation for certain three-dimen- 
sional shapes in hypersonic flow. In this “cone-expansion” 
method the flow behind the tip of a pointed body is approxi- 
mated by a Prandtl-Meyer expansion (refs. 37 and 38). 
The accuracy of this approximation is indicated by the 
comparison shown in figure 17 for the initial pressure gradient 
on an ogive of revolution. 

Tangent-cone approximation. — Newtonian impact theory 
predicts pressures depending only upon the local slope. 
This suggests approximating the pressure at each point of 
a body by that on a locally tangent cone or wedge at the 
same Mach number. For plane flow this gives equation 
(35), which yields good accuracy. For bodies of revolution, 
figure 17 gives an indication of the accuracy obtainable. 

Compression-layer approximation.- — In the upper end of 
the hypersonic range the bow shock lies close to the body 
(if the body slope is positive). This suggests making an 


approximation somewhat analogous to that of the Prandtl 
boundary-layer theory, assuming that the layer of com- 
pressed fluid between the body and shock is very thin. 

For example, assume that the shock wave lies so close 
to a circular cone that a linear variation is adequate to 
describe the stream function. Then according to equation 
(51) the stream function is given by 



(7+1 V 
2+Cy-1)ko 2 


( 1 - 0 ) 


(84) 


Requiring this to vanish at the surface gives as the ratio 
of cone angle to shock-wave angle 


5 (7+3W-2 

r 2(7+1) k 0 2 


(85) 


This result has been derived by Lees (ref. 39). At infinite 
Mach number with 7=1.405 it gives 0.916 compared with 
the true value of 0.914, and Lees shows that it is accurate 
even in the lower end of the hypersonic range. However, 
for the corresponding surface pressure coefficient, the range 
of good approximation is much smaller, as shown in 
figure 16. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 18, 1954 



APPENDIX A 


PRINCIPAL SYMBOLS 

A, B, C, D, E, T coefficients of differential equations (See 
F, G, H, I J eqs. (58) and (73).) 

B(x, y, 2 ) function defining body shape 

b(x) reduced radius (or ordinate) of axially 

symmetric (or plane) body 

b, c, d coefficients in series expansion for radius 

(or ordinate) of body (See eqs. (36) and 
(67).) 


CL 


pressure coefficient, ^ 


y 

5 


e 

k(x) 


2 P°> u « 


f, 9, h, j 


l 

m, n 

M 

P 

R{x) 

r 

Six, y, z) 

s{x) 

t 

U, V, w 
X, y, z 


0 


functions in series expansion for stream 

function (See eqs. (39) and (72).) 

initial ratio of shock-wave curvature to 

body curvature 

coefficients relating shapes 

of shock wave and body (See eq. (68).) 

free-stream Mach number 

pressure 

radius (or ordinate) of axially symmetric 
(or plane) ogive 

radius (or ordinate) in cylindrical (or plane 

Cartesian) coordinates 

function defining shape of system of shock 

waves 

radius (or ordinate) of shock wave attached 
to axially symmetric (or plane) body 
time 


V 

P 

a 

7 


* 


0 ) 

0 ) 0 , 6 > 1 , 0)2 

[] 

(") 

(~) 


velocity components in Cartesian or cylin- 
drical coordinates 

Cartesian coordinates with x in streamwise 
direction 


VM 2 -! 


( )' 
( )o 
( )» 
i ). 
( )- 


APPENDIX B 


adiabatic exponent 
semivertex angle of wedge or cone 
jO-b) 
b 

V 

conical variable, = 
x 

auxiliary hypersonic similarity parameter 
based upon local slope of shock wave, 
Mts'(x) 
r«i 

L«o 2 & <t+i) J 

bi~9b) 

9k 

2-y/2td2 (bftf fi 
37 g b 
density 

constant which is zero for plane flow, unity 
for axially symmetric flow 
thickness parameter of body; in examples, 
initial slope of shock wave 
stream function for plane or axially sym- 
metric flow (See eq. (23).) 
entropy function (See eq. (25).) 
terms in series expansion for u (See eqs. 
(40) and (70).) 

denote jump in quantity through shock 
wave 

reduced form (See eqs. (8).) 

form associated with strained coordinates. 

(See eqs. (69) and (71).) 

derivative with respect to argument 

value at tip of pointed body 

value at surface of body (or at 8, 6=b ) 

value at shock wave 

value in free stream 


CONNECTION BETWEEN HYPERSONIC AND LINEARIZED 
SUPERSONIC SIMILITUDE 


The similitude for linearized supersonic flow is now well 
understood, having been first correctly stated by Gothert in 
reference 15 (for the analogous case of linearized subsonic 
flow) . This similarity rule implies that the reduced coordi- 
nates x, y, and s of equation (8a) may again be introduced, 
and that then the reduced flow quantities, 



___1 _v_ 

T 


(Bla) 


1 w 

w= 

T 


P 

P 


M* 


vCl 0 

(£-0 


(Bib) 


depend only on the reduced coordinates and the supersonic 
similarity parameter /Jr. 1 The error in the theory and asso- 
ciated similarity rule is 0(r//3), in general. It may be em- 
phasized that here, as in all the similarity rules, the choice 
of reduced variables is by no means unique; an unlimited 
number of equivalent forms can be produced for example, by 
multiplying each reduced variable by, or adding to it, any 
constant multiple of powers of the similarity parameter. 
The particular forms adopted here were chosen to correspond 

> Here, in contract with the hypersonic case, the reduced variables ” and J> are definitely 
not 0(1) as r-»0 (for fixed AO- 
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as closely as possible to their hypersonic counterparts in 
equation (8b), so as to facilitate the following argument. 

For Mach numbers so large that M is effectively equal to /3, 
these results agree with those for the hypersonic case, 2 and 
this was pointed out by Tsien (ref. 5). However, it is more 
fruitful to reverse the argument, and observe that the hyper- 
sonic similitude, just as it stands, is entirely consistent with 
the linearized similitude. This is immediately apparent for 
the reduced velocities u, 1, and w which (as implied by the 
common notation) have identical forms in the two cases. 
They differ only in depending upon different parameters, 
but in hypersonic flow M and /3 are interchangeable to within 
the accuracy 0(r 2 ) of the theory, so that Mr can be replaced 
by dr to complete the correspondence. For the pressure 
and density, the hypersonic theory (eq. (8b)) shows that 


1 ( p 1 

V-71- 

1 'I 

tMV 

) v 

-y MV 

m 2 Vp„ 

yx 

) M 2 

(p-1) 


Again utilizing the fact that M and /3 are interchangeable in 
the hypersonic range, these can be rewritten 
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where the final forms depend upon the parameter 0t and are 
therefore (as implied by the notation) identical with their 
linearized counterparts in equation (21b). Thus the cor- 
respondence is complete. 

This means that the hypersonic small-disturbance theory, 
when properly interpreted according to the linearized super- 
sonic similitude, yields a first-order solution at all speeds 
above the transonic zone. The reduced problem of equa- 
tions (10) to (13) is solved for a given value of the parameter 
Mr, and then with Mr replaced by dr is interpreted in terms 
of physical variables according to 

[1 + t 2 u (x,y,z; dr)] 
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